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SUMMARY 


The propagation of sound in a converging-diverging duct containing a quasi- 
one -dimensional steady flow with a high subsonic throat Mach number is studied. 
The behavior of linearized acoustic theory at the throat of the duct is shown 
to be singular. This singularity implies that linearized acoustic theory is 
invalid. The explicit singular behavior is determined and is used to sketch 
the development (by the method of matched asymptotic expansions) of a non- 
linear theory for sound propagation in a sonic throat region. 


1. INTRODUCTION 


Observations of a correlation between axial Mach number and attenuation of 
sound radiated upstream from so-called sonic engine inlets have recently focus- 
ed attention on the acoustic behavior of variable-geometry ducts (refs. 1 and 
2). For. high-subsonic flows in these ducts, non-linear transonic effects be- 
come of major interest. In the linear case, a fully three-dimensional theory 
presents formidable computational difficulties, and a study of possible non- 
linear effects is, of course, even more complicated. Thus, it is natural, in 
undertaking such an effort, to restrict attention initially to a quasi-one 
dimensional model: the simplest case likely to lead to results of some prac- 

tical interest. Many earlier authors have studied linear quasi-one dimensional 
duct acoustics (see refs. 2-6, for example), but, in general, these studies 
have not been concerned with either the behavior or the validity of the lin- 
er ized solution as the axial Mach number approaches unity. 

The present paper presents some results of an ongoing analytical study of 
quasi-one dimensional acoustics in converging-diverging ducts with high-subsonic 
throat Mach numbers. The problem is inherently nonlinear, much like steady 
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transonic flow theory, but the nonlinear behavior occurs only in a narrow re- 
gion surrounding the throat section. Linearized theory yields a singular solu- 
tion in this region, and the current study is employing the method of matched 
asymptotic expansions to determine the proper solution. However, in order to 
apply any asymptotic method to correct a solution which is not uniformly valid, 
it is necessary to know in detail the singular behavior of the defective solu- 
tion. 


Thus, it is the major purpose of this paper to study the nature of the 
singularity of linearized theory at a sonic throat. To correct the defect 
using asymptotic methods requires an intricate analysis; the problem depends 
crucially on two small parameters, and the nonlinear correction to the defect 
in linearized theory involves a distinguished limit for small values of these 
parameters. The results obtained here are necessary preliminaries in this 
analysis. However, they also are of independent interest and do not appear to 
have been discussed previously. 

The analysis presented in sections 3 and 4 naturally suggests that lin- 
earization is inappropriate in a small region near the throat of the duct. The 
detailed results concerning the singular behavior of the linear solution lead 
to an appropriate stretching of the space variable and a corresponding inner 
expansion of the dependent variables which does not suffer a singularity at the 
throat. In the final section of this work the equations describing this inner 
nonlinear theory are presented, although the details of the expansion process 
are omitted for brevity. The solution of these nonlinear equations is the 
subject of current research and will appear in subsequent publications. 


2. FORMULATION AND ACOUSTIC PERTURBATION 


We consider the propagation of sound in a variable area duct carrying a 
homentropic inviscid ideal gas flow. The acoustic wavelength is assumed suffi- 
cently large, and the area variation sufficiently slow, that the field can be 
described by the equations of quasi-one dimensional gas dynamics (ref. 7): 


c p + up +pu +pu(A'/A) = 0 

St X X 

CgUf + uu x + (1/p) P x = 0 


( 2 . 1 ) 


p/p Y = constant = B 

In equations (2.1) p, p, and u are the total fluid pressure, density, and axial 
velocity, and A(x) is the duct cross sectional area. The dimensionless inde- 
pendent variables x and t are measured in units of L and L/ c s respectively , 
where L is a characteristic length associated with the area variation, and c s 
is the stagnation value of sound speed in the gas. The geometry of the problem 
is as indicated in figure 1 where the origin of x corresponds to a throat: 
A'(0)=0. 
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If the velocity and density in the basic steady flow in the duct are de- 
noted by U(x) and R(x) respectively, then from (2.1), 

UR' + RU' + RU(A' /A) = 0 , and UU' + yBR Y ~ 2 R' = 0 (2.2) 

where the energy relation in equations (2.1) has been used to eliminate the 
pressure from the system. We intend to seek solutions to the system (2.1) which 
are small perturbations about the steady values U and R, and it is convenient at 
the outset to define dimensionless variables u(x,t) and P(x,t) according to 

u(x,t) = U (x) [ 1+u (x , t ) ] , and p(x,t) = R(x) [1+p (x,t) ] (2.3) 

Substituting (2.3) into (2.1) and employing the steady relations (2.2) yields 
the system of equations on u and p in the form 

G% + M[(l+u)p +(l+p)u ] = 0 

t XX ( 2 .4) 

gSi + M(H-u)u x +(l/M)(l+p) Y_2 p x +(M , /G)[(l+u) 2 -(l+p) Y_1 ] = 0 


In equations (2.4), M(x) is the flow Mach number U(x)/c(x), c(x) is the speed 
of sound in the steady flow (c2=yBR^~^-) , and 

G(x) = (c /c) 2 = 1 + (y-l)M 2 /2 (2.5) 

s 

the latter expression following from the Bernoulli relation implied by the sec- 
ond of equations (2.2). Equations (2.4) are equivalent to those used by Cheng 
and Crocco in reference 3. 

We introduce a small dimensionless parameter 6, which measures the strength 
of the source of sound in the duct, and is assumed given from the boundary con- 
ditions associated with the system (2.1). Then u(x,t)=u(x,t;6) , p (x,t)=(x,t ; 6) 
which we assume to have expansions for 6<<1 of the form 

u = Sp(x,t) + ... , p = 6r(x,t) + ... (2.6) 


Substituting (2.6) into (2.4) and neglecting all but first-order terms we obtain 
the linearized acoustic equations: 


G 2 r 


_ + M(r +y ) = 0 
t x x 


(2.7) 


GV + Mu + (1/M)r +(M'/G)[2u-(Y-l)r] = 0 

L A X 


Equations (2.7), subject to appropriate boundary conditions, generally must 
be solved numerically because of their variable coefficients. It is the purpose 
of the present work to analyze the behavior of solutions to (2.7) in the vicin- 
ity of the throat of the duct when the throat Mach number M(0) is close to unity. 
It is well known that the system (2.7) is singular at any point where M(x)=l. 

This can be seen most simply by subtracting the two equations; the resulting 
equation has no y x term, and the coefficient of r x becomes (M 2 -l)/M, which 
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vanishes as M-KL. This can only occur at x=0 for the duct of Fig. 1. The sin- 
gularity at x=0 implies that, in general, the acoustic quantities r and y will 
be singular when the flow is sonic there. Thus, as we shall see in what follows, 
r and y generally become arbitrarily large near x=0 as M(0) approaches unity, 
thereby violating the assumptions made in deriving (2.7) that y, y , r, and r x 
all remain bounded. 

As a result of the singular behavior of the system (2.7) for high subsonic 
Mach numbers in the throat region, linearized acoustic theory is inadequate to 
describe sound propagation in the duct; we must re-formulate the perturbation 
scheme to take into account nonlinear terms in the system (2.4) which were neg- 
lected in (2.7). However, in order to make progress in this direction it is 
necessary to know precisely the nature of the singular behavior of the solutions 
y and r to (2.7). This behavior has been recognized, but never resolved, in 
previous treatments of the system (2.7) (refs. 4,5). In section 4 of this paper 
we construct an analytical general solution for the linear system (2.7) which 
displays explicitly the nature of its solutions at x=0 when M(0) is near unity. 
Before constructing this solution, however, we must discuss the behavior as 
M(0)->1 of the solutions to the steady flow equations (2.2) in some detail in 
the following section. 


3. BASIC STEADY FLOW 


As we have seen, the acoustic equations of motion are singular at x=0 when 
M(0)=1. It is useful, therefore, to introduce a parameter e=l-M(0) into our 
discussion and to consider both the basic steady flow quantities and the acous- 
tic quantities as functions of e; i.e., U(x)=U(x;e), y (x, t)= y(x, t ; e) , and so on. 
The parameter e can be considered as having been introduced through the un- 
stated boundary conditions on the steady flow. 

The elementary equations of quasi-one dimensional flow (2.2) are discussed 
in detail in numerous texts; for example, a particularly comprehensive treat- 
ment is given by Crocco (ref. 8). It is straightforward to express any of the 
fluid quantities in terms of the duct area A(x) or, equivalently, in terms of 
the Mach number M(x;e). However, the behavior of M explicitly as a function of 
x and e does not, to our knowledge, appear in the literature, and it is the pur- 
pose of the present section to determine this. 

We begin with the well-known relation implied by equation (2.2), 

M'= -MGA' / (1-M 2 )A (3.1) 

which becomes, after integration, 

cx S M s (x)[l+(y-l)M 2 (0)/2] = M s (0) [l+(y-l)M 2 (x) /2] (3.2) 

where, in equation (3.2) we have defined, 

A(x)/A(0) = a (x) , s = 2 (y-1) / (y+1) (3.3) 
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Figure 2 shows a sketch of typical integral curves of equation (3-1) assuming 
A' ’ (0)^0. We are interested in a curve such as AB in figure 2, for which M re- 
mains less than unity for all x. Since 1-M(0)=e is assumed small it is natural 
to seek an expansion of M(x;e) in the form 

M(x;e) = M q (x) + eM^(x) + ( x ) ••• (3.4) 

where Mq(0)=1. Substituting (3.4) into (3.2) and equating like powers of e we 
find that Mq(x) must satisfy 

(y+l)a s M® /2 = l+(y+l)M 0 2 /2 (3.5) 

while M^(x)=0, and 

M 2 (x) = - 2M 0 [l+(y-l)M 0 2 /2]/(y+l)(l-M 0 2 ) (3.6) 

Obviously, as can also be inferred from the integral curves of figure 2, the 
expansion (3.4) is not uniformly valid near x=0: the third term is as large as 

the first whenever 1-Mq(x) is as small as e. 

It remains to find Mq(x)=M(x; 0) in terms of x; i.e., to solve equation 
(3.5). We express a(x) as a power series 

2 

a(x) = 1 + ax + . . . 

where we assume a=A' ' (0)/2A(0)^0. Then Mq(x) can be determined after some 
algebra in the form: 

M Q (x) = 1 -((y+l)a/2) J5 |x| + ... (3.7) 

Thus, the leading term of M(x;e) behaves as a piecewise linear function of x 
near the throat so long as a^O. If a=0 we find. that Mq(x) is smooth at x=0, 
but this case will not be discussed further in this paper. 


4. SINGULARITY OF THE ACOUSTIC SOLUTION 


We shall now analyze the acoustic equations (2.7) in order to exhibit ex- 
plicitly the singular behavior of their solutions in the vicinity of the throat 
as M(0) approaches unity. Since the steady flow depends on the parameter 
e=l-M(0), the coefficients in the acoustic equations and hence the acoustic 
quantities y and r are functions of e. For e<<l we look for solutions of the 
acoustic equations in the form 

r=r Q (x,t)+er 1 (x,t)+ ... , y=y Q (x, t)+ey 1 (x, t)+ ... (4.1) 

Inserting equation (4.1) into equation (2.7) and using expansion (3.4) for the 
coefficients, we get, after neglecting higher order terms. 
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(4.2) 


G o‘ r o + V r o +M o > ‘ 0 

t XX 


G 0 y 0 + M O y 0 + + ^ M '(/ G 0^ 2y 0 ^ Y_1 ^ r (P 0 

t x x 


We recall that Mq (x) is any Mach number distribution which yields a sonic velo- 
city at the throat and a subsonic velocity throughout the remainder of the duct. 
These equations are singular at x=0 since Mq( 0)=1 *. This is most directly 
seen by subtracting the two equations. 


Analytical solutions of the system (4.2) cannot be found for arbitrary 
Mq(x) and arbitrary time dependence. However for harmonic time dependence the 
system can be reduced to a system of ordinary differential equations with a 
singular point at x=0 and no other singular points within the duct. The nature 
of this singular point will determine the singular behavior in the time harmon- 
ic acoustic quantities. Explicit analytical results concerning the exact 
nature of the singular point and the dependence on Mq(x) or Ug(x) can be found 
by use of series solution methods for linear ordinary differential equations. 

We do not present the general results of this analysis in the current paper. 
Instead, we shall illustrate the singular behavior of a general solution of the 
system (4.2) corresponding to a specific steady flow. We assume that the time 
dependence is harmonic and that the steady velocity distribution is given by a 
piecewise linear function of x: 

Ug (x) = c*(l-K|x[) |x| < (1/K) (4.3) 

where c* is the critical sound speed and K is a positive constant. This velo- 
city distribution corresponds to a reasonably shaped duct with A' (0)=0 and 
A' ' (0)/0. Equation (3.7) leads us to observe that for any duct with A' (0)=0 
and A' '(0)^0, Uq(x) will be given by equation (4.3) for x sufficiently close to 
the throat. Thus results of this section will be generally applicable to the 
throat region of many ducts of practical interest. The Mach number distribu- 
tion associated with equation (4.3) is 

Mq (x) = (1-k|x|)[(y+1)/2 - (y- 1) (1-Kj x| ) 2 /2]~^ (4.4) 

A general solution of the system (4.2) with Mq(x) given by equation (4.4) 
can be constructed by judicious use of an analytical solution found by Crocco 
and Cheng (ref. 3). In effect, they obtained a general solution to the system 
(4.2) with 


r Q (x,t) = r Q (x)exp(ifit) , y 0 (x,t) = (x)exp (ifit) (4.5) 

where = EK(2/y+l; and 

U Q (x) = c*(l+Kx) , - (1/K) <x (4.6) 

*They are also singular since Mq' has a jump discontinuity at x=0. This does 
not affect the dominant singular behavior in the solution and would not be pre- 
sent if we had chosen A' '(0)=0. 


866 



Their treatment involves introduction of the new independent variable z=(l+Kx) 2 . 
Inserting equations (4.5) and' (4. 6) into the system (4.2) and introducing z and 
o(z)=r,o(x(z)) , v(z)=pq(x(z) ) we obtain, after elimination of v(z), 

z(l-z)d 2 a/dz 2 -2[l+i3/(Y+l)]z dcr/dz - [iB(2+iB)/2(y+l)]a=0 (4.7) 

and 

v(z) = [ (y-l+iB)o-(y+l) (1-z) da/dz] / (2+iB) (4.8) 

Equation (4.7) is a hypergeometric Bquation with complex coefficients. Two 
linearly independent solutions are 

o^=F(a,b,d;l-z) and ^F(-b,-a, 2-d;l-z) (4.9) 

where 


d = 2 + 2iB/(y+l) , a+b+1 = d , ab = iB(2+iB)/2(y+l) (4.10) 

and F is the standard hypergeometric function (ref. 9). 

Since the velocity distribution used by Cheng and that in equation (4.3) 
are identical for x<0, equation (4.9) provides a general .solution to our prob- 
lem for x<0. Of course, Cheng's solution for x>0 corresponds to a supersonic 
steady flow and is not relevant to our discussion. In order to obtain a solu- 
tion when Uq(x) is given by equation (4.3) for x>0 we observe that the acoustic 
equations in this case reduce to equations (4.7) and (4.8) if 8 is replaced by 
-8. Thus we have found a general solution to the acoustic equation correspond- 
ing to U 0 (x) given by equation (4.3) for both x<0 and x>0. 

The singularity at x=0 can be found explicitly by examining the solutions 
and o 2 of equation (4.9) and the corresponding functions when B is replaced 
by -B. Clearly o± i S analytic at z=l(x=0) and the singular behavior is due to 
02> The leading term in a 2 for z near unity is given by 

o^~(l~z)^ ^F(a,b,d;0) = [cos(q£n(l-z)) + i sin(q£n(l-z)) ]/ (1-z) (4.11) 

where q=2B/ (y+1) and the - and + signs hold for x<0 and x>0, respectively (ref. 
9) . For general acoustic boundary conditions both and 0£ appear in the 
acoustic solution, and the amplitudes of the acoustic quantities will approach 
infinity as x~l when x-K). In addition their phases have an oscillatory dis- 
continuity at x=0. Figure 3 shows an example typical of the behavior of the 
acoustic quantities (pressure, in this case) for small e. The rapid rise in 
the vicinity of the throat is indicative of the developing singularity in the 
linearized acoustic quantities as e->0. For the typical case shown a pressure 
wave of magnitude unity was incident from the left on a converging-diverging 
section situated in an otherwise uniform duct. 

Since any duct with a tnroat at which A'*(0)f0 will have a locally (near 
x=0) piecewise linear steady velocity distribution, the acoustic quantities in 
such a duct will have the singular behavior given in equation (4.11). In this 
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circumstance the linear theory governing expansion (2.6) fails for any nonzero 
<5, no matter how small. Nonlinear effects become appreciable in the throat 
region. In the next section, an approach is outlined which applies the method 
of matched asymptotic expansions to study the nonlinear effect in the region 
near the throat. 


5. NONLINEAR PERTURBATION EQUATIONS 


In this section we set forth in summary the theory describing sound propa- 
gation near the throat as M(0)-KL. We regard the expansion indicated by equa- 
tions (2.3) and (2.6) as an outer expansion valid as <5->0 for fixed values of x, 
t , and e ; i . e . , 


u(x,t;e,<5) = U(x;e)[l + <5p(x,t;e) + ... ] 
p(x,t;e,6) = R(x;e)[l + 6r(x,t;e) + ... ] 


(5.1) 


From the details of the previous section we know that in general p and r become 
arbitrarily large in the limit as x and e approach zero, being expected to grow 
as x~l. 

__ Therefore we introduce an inner variable X=[ (y+l)/2] 35 (x/e) and assume that 
u and p behave asymptotically as e-K) with X, t fixed as: 

u = u 1 (X,t;e) = U i (X;s)[l+ey 1 (X,t) + . . . ] 

. , . (5.2) 

p = p (X,t;e) = R (X;e)[l+er x (X,t) + ...] 

where e is assumed to be a function of 8 which vanishes as <5~K) and is to be 
determined by asymptotic matching of (5.2) with (5.1). In addition we expand 
the steady flow quantities in the form: 

U 1 (X;e)=uJ(X)+eU^(X)+... , R 1 (X;e)=Rj(X)+eR^ (X)+. . . (5.3) 

Equations (5.3) could be combined directly with (5.2) as one inner expansion, 
but we find that it simplifies the considerable algebra involved to retain the 
dimensionless perturbations p and r and to multiply the separate expansions as 
indicated in equation (5.2). 

The steady flow quantities in equation (5.3) are known if the correspond- 
ing expansion of the Mach number is known. We assume an expansion for M of the 
same form as (5.3), substitute this expansion into equation (3.1), transform to 
the inner variable X, and solve the successive differential equations which 
result. After matching with the outer expansion (3.4) we find 

M(x;e) = M 1 (x;e) = 1 - e (1+aX 2 )^ + ... (5.4) 

where a=A' ' (0) /2A(0) as before. Using (5.4) we find the coefficients in the 
expansions (5.3) by use of the steady flow relations between M and U or R. 
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Next we carry out the same process on system (2. A) using u=ep 1 + 
p=er 1 +. . . , and substituting equation (5.4) in the coefficients. This yields 
the system of inner equations satisfied by p 1 and r 1 in the form: 


' 2 35 

n t + U+n-(i+aX z ) ] tt x - 
?x = 0 


- (c+n) » 0 

(1+aX 2 ) 


(5.5) 


where in equations (5.5) we have defined 

C = (3-y) (y i +r i )/4 and n = (y+1) (p^r 1 ) /4 


Equations (5.5) are the nonlinear equations which, to first order in e, 
govern sound propagation through a throat as the Mach number there approaches 
unity. The lengthy details of solving the system will not be presented here. 
However, certain important conclusions can be made at this stage. The quanti- 
ties q and £ are related to the Riemann invariants of system (2.4), n represent- 
ing the upstream and £ the downstream propagating portions of the solution to 
(2.4). Considerations of asymptotic matching between expansions (5.1) and (5.2) 
lead to the conclusion that, to first order in e, ? actually vanishes. Thus, 
as is often argued from physical considerations, the lowest order nonlinear 
effect of the sonic throat is on the upstream propagating waves alone. 


A final observation which we make here is that matching considerations 
indicate that, in the distinguished limit implied by the inner expansion (5.2), 
e is to be taken equal to 6^. Hence we conclude that, given an acoustic source 
strength 6, nonlinear effects on sound arise for throat Mach numbers as far 
away from unity as 8% . This would explain why marked sonic throat effects are 
observed experimentally for throat Mach numbers as low as 0.75-0.8. 
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